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Flow about the NASA Trapezoidal Wing is computed with several turbulence models by using grids from 
the first High Lift Prediction Workshop in an effort to advance understanding of computational fluid dynamics 
modeling for this type of flowfield. Transition is accounted for in many of the computations. In particular, a 
recently-developed 4-equation transition model is utilized and works well overall. Accounting for transition 
tends to increase lift and decrease moment, which improves the agreement with experiment. Upper surface flap 
separation is reduced, and agreement with experimental surface pressures and velocity profiles is improved. 
The predicted shape of wakes from upstream elements is strongly influenced by grid resolution in regions above 
the main and flap elements. Turbulence model enhancements to account for rotation and curvature have the 
general effect of increasing lift and improving the resolution of the wing tip vortex as it convects downstream. 
However, none of the models improve the prediction of surface pressures near the wing tip, where more grid 
resolution is needed. 


Nomenclature 


a speed of sound 

C f surface skin friction coefficient, 2 r w / (pU^ e f ) 

C p surface pressure coefficient, 2 (p — p re f)/ {pU^ e f) 

Cd drag coefficient 

C l lift coefficient 

Cm moment coefficient about y- axis 

E error compared to experiment 

e N transition prediction method based on linear stability theory 

FT fully turbulent 

i t turbulence index, \dv / dn\ / (nu T ) 

k turbulence kinetic energy (TKE) 

L reference length 

M freestream Mach number, U re f / a re f 

n normal direction to wall 

p pressure 

Re freestream Reynolds number, U re fL/ u re f 

Reo momentum- thickness Reynolds number, XJQ/v 

Tu turbulence intensity (percent), 100^2/3(fc/L^? e ^) 

V velocity magnitude 

Ui (u,v,w) Cartesian velocity component 
u T wall friction velocity, yjr w / p 

Xj (x,y,z) Cartesian coordinate 
a angle of attack, deg. 

d partial derivative 

6 momentum thickness 

7 intermittency parameter 
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K 

Von Karman law-of-the-wall constant, 0.41 

Pt 

turbulent eddy viscosity 

P 

density 

V 

kinematic viscosity 

V 

turbulent variable from Spalart-Allmaras model 

T ij 

turbulent shear stress tensor 

T~w 

wall shear stress 

n 

vorticity magnitude 

Subscript 


max 

maximum 

ref 

reference conditions (free stream) 

X 

based on ^-distance 


I. Introduction 

The first AIAA High Lift Prediction Workshop (HiLiftPW-1) established a baseline for current state-of-the-art 
computational fluid dynamics (CFD) code capabilities on high-lift configurations . 1,2 At the workshop, a different 
research groups using a variety of codes and methods studied the 3-element NASA Trapezoidal Wing configuration . 3,4 
Particular assessments were made for fully turbulent flow with two different flap settings. The importance of transition 
was known, but the first workshop’s cases were limited to fully turbulent flow to attain more useful comparisons 
between a large number of codes. As a whole, the collective CFD results on a clean configuration (no slat or flap 
support brackets modeled) tended to under-predict the lift and the magnitude of the pitching moment compared to 
experiment. However, the trend with grid refinement was generally in the right direction (i.e., finer grids produced 
higher lift). Including slat and flap support brackets in the CFD geometry tended to lower the predicted lift. Agreement 
(both between different CFD results as well as between CFD and experiment) deteriorated at higher angles of attack 
approaching stall. The outboard region of the wing was the most problematic area for CFD, with the least consistency 
among CFD codes, and the largest deviations from experiment. 

There were many publications of CFD results after the workshop . 5-18 Most of these studies involved a few well- 
known “standard” turbulence models such as Spalart-Allmaras (SA ) 19 and Menter’s shear-stress transport (SST ) 20 
applied in fully turbulent mode. At the workshop, SA (fully turbulent) was heavily favored because it tended to yield 
closer agreement to experiment for the Trapezoidal Wing. SST tended to separate more readily than SA and con- 
sequently produced lower lift. However, it was noted at the workshop that transition effects could alter this trend. 
Transition tended to reduce trailing edge flap separation and significantly improve pitching moment predictions. Sub- 
sequently, the papers of Steed , 6 Fares and Nolting , 11 and Eliasson et al . 9, 17 explored the important effects of including 
transition in the CFD simulations. Steed used the 4-equation 7 -Reg SST transition/turbulence model of Langtry and 
Menter . 21 This model includes the solution of two additional transport equations for an intermittency parameter and 
a local transition onset momentum- thickness Reynolds number. The effective intermittency regulates the production 
of turbulence kinetic energy (TKE) in the SST model. Fares and Nolting fixed laminar regions in a lattice Boltz- 
mann solver based on measurements described by McGinley et al . 22 Eliasson et al. used a database e N method with 
the envelope approach to impose transition regions in conjunction with the SA model. All groups demonstrated im- 
proved results when transition regions were included. The importance of including transition for multi-element airfoil 
computations was similarly shown to be very important for two-dimensional flows over a decade ago . 23 

In the current validation study, some of the best practices from the earlier studies mentioned above are applied to 
the Trapezoidal Wing in an effort to advance our understanding of high lift flowfields. Many cases are included, using 
different codes, grids, and turbulence models. Through particular cases in the matrix of runs, an attempt is made to 
demonstrate both isolated and combined effects. The study includes the 7 -Reg SST transition/turbulence model of 
Langtry and Menter , 21 and we also investigate setting transition locations in an a priori fashion. Both structured- and 
unstructured-grid solvers are employed, and results are examined both with and without support brackets included. 

The influences of rotation and curvature corrections to turbulence models are also investigated. These could be 
influential for regions of the flow such as near the wing tip where a strong vortex forms and subsequently convects 
downstream. Turbulence model variations tested include the three rotation and curvature-corrected turbulence models 
SA-R , 24 SA-RC 25 and SST-RC . 26,27 In particular, the ability of these models to resolve the Trap Wing’s tip vortex are 

a http://hiliftpw.larc.nasa.gov; cited 5/20/2012. 
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assessed compared to the baseline SA and SST models. 


II. Methods and Models 

Two CFD codes are employed in this study. CFL3D is a cell-centered, structured-grid, upwind-biased Reynolds- 
averaged Navier-Stokes (RANS) code , 28 and FUN3D is a node-centered, unstructured-grid, upwind-biased RANS 
code . 29 These codes were also used in an earlier study associated with HiLiftPW-1 . 12 Transition is accounted for in 
two different ways: ( 1 ) by zeroing out turbulence production terms in laminar regions on each of the elements by using 
regions specified by Eliasson et al . 17 (also provided on the HiLiftPW website), and (2) by employing the 4-equation 
7 -Ree SST transition/turbulence model of Langtry and Menter . 21 

In CFL3D and FUN3D, the default implementations of the production terms for the SST and 7 -Ree SST models 
make use of the “vorticity approximation,” introduced by Menter 30 and referred to elsewhere 0 as “SST-V.” In this 
approximation, the production term Tijdui/dxj is replaced with /i t £ 2 2 , where fi t is the eddy viscosity and £2 is the 
magnitude of vorticity. This approximation has been shown to be reasonable for many aerodynamic flows, and it 
avoids over-production of TKE through strong shocks and near stagnation points that occurs with the standard strain 
production form. A separate run of the 7 -Ree SST model was made with CFL3D using the standard strain production 
form, and it made a dramatic difference (4% decrease in lift, 5% decrease in drag, and 7% increase in moment) 
compared to the use of the vorticity production form. With the standard strain production form, the tendency of 
overproduction of TKE near stagnation points adversely influences the 4-equation model’s attempt to laminarize near 
the leading edges of the main and flap elements. Note that there are other production term limiters (Langtry 31 ), 
modified production term forms (Kato and Launder 32 ), and time scale limiters (Durbin 33 ) that also improve stagnation 
region behavior. All results shown below use the vorticity approximation (“SST-V”). 

Several rotation and curvature corrections to the SA and SST models were investigated. The SA-R model 24 reduces 
the eddy viscosity in regions where vorticity exceeds strain rate, such as in vortex core regions where pure rotation 
should not produce turbulence. The SA-RC model 25 accounts for both rotation and curvature effects in the SA model 
through a function of the Lagrangian derivative of the strain rate tensor. The SST-RC model 26,27 sensitizes the SST 
model to curvature based on an approximation to the gradient Richardson number. Note that for SA-RC, the constant 
c r 3 was taken as 0.6 (based on an earlier study 34 ), and for SST-RC the constant C rc = 1.4 was used from Mani et al . 27 

As one of several validation exercise for the recent implementation of the 7 -Ree SST model in NASA codes, 
the flat plate was computed using various freestream turbulence intensity and eddy viscosity levels. The cases were 
taken from Langtry and Menter . 21 Here, the codes CFL3D, FUN3D, and OVERFLOW 35 were computed on the same 
hexahedral grid with plate length of 4.367L, M = 0.2, and Re/L = lx 10 6 . Skin friction results are shown in Fig. 1 
on page 1 1 . All three codes yielded similar results to each other, and agreed reasonably well with experiment. 

III. Results 

For the NASA Trapezoidal Wing, computations were performed for configuration 1 (slat at 30°, flap at 25°). 
The Re based on mean aerodynamic chord ( L = 39.634 inches) was 4.2 million and M = 0.2. A range of angles 
of attack between 6 ° and 37° were tested, with a focus on the a = 13° and a = 28° cases. Table 1 shows a 
summary of computations performed for this study. In the 7 -Ree SST runs, the far-field turbulence intensity was set 
to Tu = 0.075% (from NASA Langley 14- by 22-foot wind tunnel conditions 36 ), and far-field eddy viscosity was 
set to fi t / fi re f = 0.2. A higher value of fi t /fi re f = 10 was also tried for one case, and it made little difference. 
So-called “sustaining terms ” 37 were not employed for the current study, due to uncertainty in their effect on the 7 -Ree 
SST model, which had been calibrated by Langtry and Menter without them. The term “FT” in the table indicates 
fully turbulent, and the term “trans set” indicates that transition locations were fixed a priori, as described in Section 
II. Note that some of the listed computations were actually a part of the earlier study by Park et al ., 12 but they are 
included here for comparison purposes. 

To provide some representative results in tabular form, a summary of selected forces and moments for various 
run combinations of code, grid, turbulence model, and transition treatment at a = 13° and 28° is given in Table 
2. Note that it proved to be very difficult to run at higher angles of attack such as a = 28° with transition fixed 
a priori. The flowfield tended to separate (stall) early and resist convergence. Fixing transition locations was also 
tedious and required knowing locations prior to running, which is not generally feasible unless coupled to a method 
such as e N . Furthermore, the regions where laminar flow was specified sometimes appeared to retain too-high levels 

b http://turbmodels.larc.nasa.gov; cited 5/20/2012. 
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of turbulence (because of turbulence advection from upstream). Various options (such as SA-R) could mitigate the 
degree of upstream influence, but the process in general was not very straightforward. 

Three different grid systems, available on the HiLiftPW website, were employed in this study. All computations 
were performed on “medium” grid levels from the workshop grid refinement study. Grid study results for the UH6 and 
SX1 grids can be found in Park et al. 12 The grids denoted as UH6 (without brackets) and UH16 (both with and without 
brackets) were unstructured mixed-element grids. SX1 (without brackets) was a structured multi-zone grid, and UX9 
was the same grid written in unstructured format (hexahedra). For the medium levels, the mixed-element UH6 grid 
(no brackets) had 11 million grid points, UH16 (no brackets) had 28 million, UH16 (brackets) had 32.5 million, and 
SX1/UX9 (no brackets) had 52 million. Additional details about these grids can be found in other references. 2,5, 13 At 
HiLiftPW- 1, including the brackets was found to have a noticeable influence on results. For example, at a = 28° the 
lift coefficient decreased by approximately 0.06-0.08 when brackets were modeled. 

Figure 2 on page 1 1 shows two representative views of the underside of the configuration, one without and one 
with brackets. Figure 2(a) is the structured SX1 grid, and (b) is the unstructured UH16 grid. The surface grid points 
on the main element are shown in red. Figure 3 on page 12 shows the wing upper surface, along with locations where 
pressure taps on the model were located and where 7-hole-probe velocity data were taken. In the current study, surface 
pressures are compared only at 3 stations on the flap, as depicted in the figure by orange. 

Table 1: Summary of computations 


Code 

Grid 

Brackets? 

Turb Model 

Focus 

FUN3D 

UH6 

No 

SA (FT) 

full polar 

FUN3D 

UH16 

No 

SA (FT) 

full polar 

FUN3D 

UH16 

Yes 

SA (FT) 

full polar 

FUN3D 

UH16 

Yes 

SA-RC (FT) 

limited a 

FUN3D 

UH16 

Yes 

SA (trans set) 

limited a 

FUN3D 

UH16 

Yes 

SA-R (trans set) 

limited a 

FUN3D 

UH16 

Yes 

SA-RC (trans set) 

limited a 

FUN3D 

UH16 

Yes 

SST-RC (trans set) 

limited a 

FUN3D 

UH16 

Yes 

7 -Re e SST 

full polar 

FUN3D 

SX1/UX9 

No 

SA (FT) 

limited a 

CFL3D 

SX1/UX9 

No 

SA (FT) 

full polar 

CFL3D 

SX1/UX9 

No 

SST (FT) 

full polar 

CFL3D 

SX1/UX9 

No 

SA-R (FT) 

limited a 

CFL3D 

SX1/UX9 

No 

SA-RC (FT) 

limited a 

CFL3D 

SX1/UX9 

No 

7 -Ree SST 

full polar 

CFL3D 

SX1/UX9 

No 

7 -Re e SST-RC 

full polar 


A. Forces and Moments 

Figure 4 on page 13 shows a summary of the computed lift, drag, and pitching moment coefficients for several different 
representative cases. CFL3D was only run without brackets. As shown in the left-half column (Figs. 4(a), (c), and (e)), 
for CFL3D the SA model run fully turbulent (“FT”) produced very reasonable results compared to experiment, slightly 
low in lift over most of the lift curve. The SST model run fully turbulent yielded significantly lower lift — especially 
at high angles of attack — and a higher moment overall. As noted at the workshop, SST tended to produce a greater 
amount of separated flow than SA on the flap for many workshop participants. On the other hand, the 4-equation 
7 -Reg SST model increased the lift and decreased the moment, in excellent agreement with the experiment overall. 
Below we will explore some detailed results from the 7 -Reg SST model in order to try to understand why this was 
the case. However, this particular model also resulted in Cl, max that was too high on the SX1 grid. The reasons for 
this are not known, although including the brackets may be necessary to accurately capture some of the flow physics 
mechanisms responsible for stall. 
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Table 2: Summary of selected forces and moments 


Code 

Grid 

Brackets? 

Turb Model 

C L 

Cd 

Cm 

%Ecl 

%Ecd 

%Ecm 

a = 13° 










Experiment 




2.0474 

0.3332 

-0.5032 




FUN3D 

UH6 

No 

SA (FT) 

1.9922 

0.3248 

-0.4707 

-2.7 

-2.5 

-6.5 

FUN3D 

UH16 

No 

SA (FT) 

1.9939 

0.3254 

-0.4714 

-2.6 

-2.3 

-6.3 

FUN3D 

UH16 

Yes 

SA (FT) 

1.9801 

0.3261 

-0.4676 

-3.3 

-2.1 

-7.1 

FUN3D 

UH16 

Yes 

SA-RC (FT) 

1.9886 

0.3276 

-0.4698 

-2.9 

-1.7 

-6.6 

FUN3D 

UH16 

Yes 

SA (trans set) 

2.0132 

0.3317 

-0.4812 

-1.7 

-0.5 

-4.4 

FUN3D 

UH16 

Yes 

SA-R (trans set) 

2.0226 

0.3350 

-0.4857 

-1.2 

0.5 

-3.5 

FUN3D 

UH16 

Yes 

SA-RC (trans set) 

2.0454 

0.3354 

-0.4938 

-0.1 

0.7 

-1.9 

FUN3D 

UH16 

Yes 

SST-RC (trans set) 

2.0271 

0.3289 

-0.4816 

-1.0 

-1.3 

-4.3 

FUN3D 

UH16 

Yes 

7 -Ree SST 

2.0122 

0.3309 

-0.4797 

-1.7 

-0.7 

-4.7 

FUN3D 

SX1/UX9 

No 

SA (FT) 

1.9974 

0.3232 

-0.4771 

-2.4 

-3.0 

-5.2 

CFF3D 

SX1/UX9 

No 

SA (FT) 

2.0165 

0.3254 

-0.4831 

-1.5 

-2.3 

-4.0 

CFF3D 

SX1/UX9 

No 

SST (FT) 

1.9848 

0.3180 

-0.4664 

-3.1 

-4.6 

-7.3 

CFF3D 

SX1/UX9 

No 

SA-R (FT) 

2.0295 

0.3274 

-0.4862 

-0.9 

-1.7 

-3.4 

CFF3D 

SX1/UX9 

No 

SA-RC (FT) 

2.0434 

0.3305 

-0.4923 

-0.2 

-0.8 

-2.2 

CFF3D 

SX1/UX9 

No 

7 -Re e SST 

2.0581 

0.3312 

-0.4996 

0.5 

-0.6 

-0.7 

CFF3D 

SX1/UX9 

No 

7 -Re e SST-RC 

2.0732 

0.3345 

-0.5072 

1.3 

0.4 

0.8 

P 

II 

to 

00 

o 










Experiment 




2.8952 

0.6776 

-0.4559 




FUN3D 

UH6 

No 

SA (FT) 

2.8478 

0.6692 

-0.4257 

-1.6 

-1.2 

-6.6 

FUN3D 

UH16 

No 

SA (FT) 

2.8470 

0.6650 

-0.4250 

-1.7 

-1.9 

-6.8 

FUN3D 

UH16 

Yes 

SA (FT) 

2.7950 

0.6537 

-0.4133 

-3.5 

-3.5 

-9.3 

FUN3D 

UH16 

Yes 

SA-RC (FT) 

2.7936 

0.6581 

-0.4053 

-3.5 

-2.9 

-11.1 

FUN3D 

UH16 

Yes 

7 -Reg SST 

2.8674 

0.6746 

-0.4390 

-1.0 

-0.4 

-3.7 

FUN3D 

SX1/UX9 

No 

SA (FT) 

2.8594 

0.6650 

-0.4288 

-1.2 

-1.9 

-5.9 

CFE3D 

SX1/UX9 

No 

SA (FT) 

2.8457 

0.6613 

-0.4319 

-1.7 

-2.4 

-5.3 

CFE3D 

SX1/UX9 

No 

SST (FT) 

2.8073 

0.6476 

-0.4186 

-3.0 

-4.4 

-8.2 

CFE3D 

SX1/UX9 

No 

SA-R (FT) 

2.8573 

0.6632 

-0.4358 

-1.3 

-2.1 

-4.4 

CFE3D 

SX1/UX9 

No 

SA-RC (FT) 

2.8915 

0.6754 

-0.4472 

-0.1 

-0.3 

-1.9 

CFE3D 

SX1/UX9 

No 

7 -Reg SST 

2.9010 

0.6743 

-0.4519 

0.2 

-0.5 

-0.9 

CFE3D 

SX1/UX9 

No 

7 -Re e SST-RC 

2.9038 

0.6756 

-0.4521 

0.3 

-0.3 

-0.8 
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In the right-hand column (Figs. 4(b), (d), and (f)), FUN3D results are shown for the UH16 medium grids both 
without (“no Brkts”) and with (“Brkts”) brackets. FUN3D’s SA results on the configuration without brackets were 
very similar to those of CFL3D. Although not shown for clarity, the FUN3D force and moment results for the UH 6 
grid were very similar to the UH16 fully turbulent SA results through stall. Using SA fully turbulent with brackets 
resulted in Cl, max ~ 2.82, around 6 % too low. Use of the 4-equation 7 -Ree SST model resulted in lift levels on the 
grid with brackets very close to experiment (similar to SA fully turbulent on the grid without brackets). However, at 
a = 34°, the 7 -Ree SST solution stalled (too early) on the grid with brackets and did not converge. 

Thus, when computing on a grid that includes brackets, CFD results tended to more closely approach experimental 
levels over most of the lift curve when transition was included. As understood from these results (in concert with 
earlier workshop studies mentioned in the introduction), there are competing mechanisms at play when attempting 
to predict global forces and moments for the NASA Trapezoidal Wing. In terms of lift: including brackets tends to 
decrease lift, refining the grid tends to increase it, and including transition tends to increase it. Turbulence model 
also has a clear effect, with SA tending to produce less separation than SST, all other things being equal. Additional 
turbulence model influences will be explored below. 

B. Surface Pressures 

Figure 5 on page 14 and Fig. 6 on page 15 show surface pressure coefficients at selected stations for a = 13° and 
28° by using CFL3D and FUN3D, respectively. At the upper surface flap forward station, there was a trend of lower 
pressures (more negative C p ) for the 7 -Ree SST model, followed by the SA model, and higher pressures for the SST 
model. Generally, the results from the 7 -Reg SST model agreed best with experiment along this station, and fully 
turbulent SST the worst. The FUN3D results on the grid with brackets captured the “scalloping” effect of the brackets 
on the pressures, in better agreement with experiment. All models performed poorly near the wing tip region, outboard 
of y « —80 inches. This was also the case at the first workshop, and was most likely due to inadequate resolution of 
the wing tip vortex. Either insufficient grid density or deficiencies in turbulence modeling (or a combination of both) 
was probably to blame. A subsequent study by Lee and Pulliam 18 demonstrated the important influence of grid density 
in the wing tip region. 

At the 50% span station on the flap, the 7 -Ree SST model yielded results closer to experiment than fully turbulent 
SA or SST at a = 13°. The tendency of this transition model to yield lower pressures over the flap upper surface 
can clearly be seen. At a = 28°, none of the models agreed well with experiment. However, the FUN3D results 
were somewhat closer to experiment than those of CFL3D, probably because of the bracket influence at this particular 
station . 1 

Finally, at the 85% span station on the flap, similar trends were seen for the three turbulence models. However, 
in some cases the 7 -Ree SST model agreed best with experiment, while in others the fully turbulent SA model was 
closer. The brackets also had an influence at this span station: including them improved the comparison with lower 
surface flap pressures at both a = 13° and 28°. 

C. Off-Body Velocity Profiles at a = 28 ° 

Recently, off-body velocity profiles (measured using 7-hole probes) at a limited set of locations in the NASA Langley 
14- by 22-foot wind tunnel were made available in Hannon et al . 4 One set of these experimental data were taken at 
a = 28°. Velocity magnitude results on the SX1/UX9 grid are shown in Fig. 7 on page 16, using SA (fully turbulent), 
SST (fully turbulent), and 7 -Reg SST at four locations normal to the wing and flap upper surfaces. These surface 
locations are near: x = 17.5, y = —13.6; x = 64.2, y = —71.5; x = 69.7, y = —71.9; x = 75.6, y = —72.4 
inches. Unfortunately, the experimental data probe did not traverse close enough to the surface over the main element 
(Figs. 7(a) and (b)), so it missed the boundary layer and most of the slat wake. CFL3D and FUN3D were both run 
with the SA model (fully turbulent) on this same grid, and their results were very similar at each station. CFD results 
at the aft flap station (in the region where surface streamlines indicated separated flow) were all poor, but at the other 
three stations the CFD results agreed fairly well with experiment, and the 7 -Ree SST model performed generally best. 
It yielded particularly good results for the wakes at the forward flap station. Experimental data were also available 
for the respective x, y, and 2 components of velocity at three of the stations, but plotting these did not reveal any 
additional information helpful for turbulence model validation. 

Figure 8 on page 17 compares several FUN3D results on the UH16 grid. Fully turbulent SA model results, on grids 
both with and without brackets, were very similar to each other. The 7 -Ree SST model results were generally better 
compared to experiment. However, all results on the UH16 grid were much worse compared to experiment (wakes 
were more diffused) than results on the SX1/UX9 grid (see Fig. 7). By inspecting cuts of the respective grids near 
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the outboard span stations (Fig. 9 on page 18), it is clear that the SX1/UX9 grid provided much greater wall-normal 
grid density than UH16 in the regions above the wing where the wakes from upstream elements pass. This greater 
grid density allowed more accurate wake resolution. Providing grids with pre-specified wake clustering is currently a 
challenge for many unstructured grid generation codes. 

Experimental off-body velocity and turbulence profile measurements are extremely useful for CFD validation 
studies such as this, because they can highlight specific deficiencies in grid resolution and turbulence modeling. As 
seen in the current velocity results, inadequate grid resolution in wake regions can produce diffuse wake structures that 
do not resemble the true physics. It is unclear what impact the diffused wakes have on global quantities of interest, 
because so many other factors also come into play. However, Smith 38 made the observation that when wakes and 
boundary layers merge too quickly, early separation will set in. This implies that under-resolution of the wakes in a 
CFD computation could influence stall prediction. 

D. Transition Effects 

Although the 7 -Rcq SST model yielded very good results for this configuration, the model did not converge well 
as implemented in both CFL3D and FUN3D. In particular, the residual of the 7 turbulence equation tended to stall, 
typically not decreasing beyond a certain level (depending on the case). An example using CFL3D at a = 13° is 
shown in Fig. 10 on page 18. Here, the mean flow density residual (red line) was converging, and lift coefficient 
(blue line) was settling to nearly a constant level. However, the residual for the 7 turbulence equation (green line) was 
oscillatory near 10 -7 , and not dropping. This represents residual convergence of only about one order of magnitude 
from this equation’s initial (maximum) level. Nonetheless, in general we have observed that the patterns of laminar 
and turbulent flow generally “converge” in the sense that they achieve a fixed pattern then do not change noticeably. 
Naturally, poor convergence of turbulence equation residuals is a source of concern and potential uncertainty because 
an adequately converged mean flow final result is not assured. More work needs to be done to investigate the numerical 
behavior of the 7 -Ree SST model. 

The difference between the flowfields in the flap region using SST and the 7 -Ree SST models is exemplified in 
Fig. 11 on page 19. The fully turbulent SST model yielded lower velocities over the crest of the flap than the 7 -Ree 
SST model (Figs. 11(a) and (b) at a representative span station), and SST produced a larger region of separated flow 
on the aft half of the flap. The nondimensional eddy viscosity levels for the two models are shown in close-ups over 
the front part of the flap at the same representative span station in Figs. 1 1(c) and (d). Generally, levels of fi t / fjb re f > 1 
are considered turbulent. The 7 -Ree SST model produced “laminar” levels of eddy viscosity (seen as thin dark blue 
regions) over about the first 40% of the flap upper surface, as well as over the entire lower surface at this ^/-station. 
According to Smith , 38 extensive laminar flow helps to delay separation because it reduces Ree. Further, thin boundary 
layers can withstand stronger adverse gradients than thick ones. Upstream transition effect on downstream wake shape 
may also play a role . 23,39 The upper surface streamlines for the two models are shown in Figs. 1 1(e) and (f). The 7 -Ree 
SST model yielded significantly less separated flow on the flap than fully turbulent SST. 

Next the predicted transition locations from the 7 -Ree SST model are compared with those of Eliasson et al . 17 To 
visualize the laminar and turbulent regions on the Trapezoidal Wing, an approximated turbulence index i t was used, 
based on the one described in Spalart and Allmaras . 19 This index was designed to yield a value close to zero in a 
laminar region and close to one in a turbulent region for the SA model. Although not designed with other turbulence 
models in mind (which may have different asymptotic behavior near the wall), translating any model’s eddy viscosity 
to an equivalent v and making use of the Spalart- Allmaras formula has been found to be valuable for visualizing 
regions where there are clear demarcations between laminar and turbulent flow. 

Figure 12 on page 20 shows i t contours compared with locations predicted by Eliasson et al. using the e N method, 
with N = 10. Eliasson et al. found this latter method (in combination with SA) to yield reasonably good correlations 
with the experimental data of McGinley et al . 22 as well as produce good overall force and moment results in compu- 
tations compared to experiment. On the upper surfaces (Figs. 12(a), (c), and (e)), there was generally good correlation 
between the methods. On the lower surfaces (Figs. 12(b), (d), and (f)), the correlation between the methods was fairly 
good on the flap (where mostly laminar flow was indicated), but was not as good on the main element. Here, the 7 -Ree 
SST model indicated either transitional or turbulent values in very unusual patterns, whereas the e N method indicated 
more significant regions of laminar flow. In contrast, McGinley et al. noted primarily turbulent flow near the mid-chord 
region of the main element lower surface in the experiment. As mentioned above, the current approximate turbulence 
index is only helpful for determining regions that are very definitively laminar or turbulent. On the lower surface of 
the main element the 7 -Ree SST model produced results that are indeterminate with this metric. Although not shown, 
i t contours for the 7 -Reg SST model in FUN3D on the UH16 grid with brackets looked similar to those in Fig. 12, 
although the presence of the brackets caused some non-uniformity in the index values on the adjacent downstream 
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element. 

In Fig. 13 on page 21, examples of the combined effects of brackets and transition on Cl at a = 13° are shown. 
With no brackets and SA (fully turbulent), the lift was computed about 2.6% too low. Adding brackets lowered lift 
(making it 3.3% low). However, including transition on the grid with brackets — either through fixing transition a priori 
or using the 7 -Ree SST model — brought the lift back up (making it 1.7% low, near the lower experimental uncertainty 
bounds). 

E. Rotation and Curvature Correction Effects 

Looking beyond the trends caused by including transition effects, we now examine the effects of including various 
rotation and curvature corrections in the turbulence models. These corrections are expected to better account for 
areas of the flow where either pure rotation should not produce additional turbulence, or where streamline curvature 
near solid walls should influence turbulence (reduction near convex curvature, enhancement near concave curvature). 
Figure 14 on page 21 shows nondimensional vorticity contours for a = 13°, for two representative cases without and 
with brackets. These plots show flowfield features for which rotation and curvature corrections may be important. In 
particular, the wing tip sheds a strong vortex that convects downstream. When brackets are included, there are other 
smaller vortical features shed inboard as well. Including rotation and curvature corrections had a significant influence 
on the strength of the shed vortices. For example, at the aft-most station shown ( x = 90 inches), peak tip vortex 
strength was increased (compared to baseline SA or SST, which were similar) by 15% with SA-R and by 21% with 
SA-RC. The 7 - Ree SST model increased tip vortex strength by 10%, and 7 -Ree SST-RC increased tip vortex strength 
by 28% over the baseline models. 

Figure 15 on page 22 summarizes the effects of the various rotation/curvature corrections. In Figs. 15(a) and (b), 
lift values at the two angles of attack of a = 13° and 28° are shown from CFL3D with different turbulence model 
variants. Use of SA-R, SA-RC, or SST-RC had the influence of increasing Cl (trends were similar for FUN3D). The 
SA-RC model had a larger influence than SA-R. The SST-RC correction was run in conjunction with the transition 
model Ci-Reg SST-RC), and yielded higher lift than 7 -Reg SST alone. Although not shown in figures, effects on drag 
and moment can be seen in Table 2. Including transition and rotation/curvature corrections both had the tendency to 
increase drag and decrease moment, usually improving correlation with experiment. 

Finally, we address the question of whether the turbulence model rotation and curvature corrections improved 
surface pressure predictions in the problematic wing tip region. Figures 15(c) and (d) show span wise surface pressure 
coefficient comparisons at the flap forward span station. Two angles of attack ( a = 13° and 28°) are shown for CFL3D 
on the SX1 grid. Generally, the rotation and curvature corrections tended to decrease the upper surface flap pressures, 
but none of the models improved the tip region predictions. As implied in Lee and Pulliam , 18 significantly greater grid 
resolution is probably required in the wing tip region. Only with adequate grid resolution will it be possible to conduct 
meaningful studies of the influence of turbulence models on the flow physics in this area of the flow. Figure 15(e) 
shows results for FUN3D on the UH16 grid with brackets. Here, for clarity only SA and SA-RC results are shown 
(both with transition set a priori). Again, the rotation/curvature correction yielded lower upper surface pressures along 
most of this station on the flap, in better agreement with experiment out to approximately y = — 65 inches, but did not 
improve the predictions near the wing tip. In Fig. 15(f) the differences in predicted C p between SA-RC and SA are 
shown in a contour plot. Most of the largest differences occurred near the wing tip. On the flap, the C p levels from 
SA-RC were mostly lower over the front part and mostly higher near the trailing edge. 

IV. Conclusions 

The first High Lift Prediction Workshop, which used the NASA Trapezoidal Wing experimental data, brought 
several CFD issues to light. Among these were the important influences of grid resolution, brackets, and transition. The 
current paper’s goal was to make use of many of the lessons learned from the workshop and subsequent independent 
studies, and attempt to further advance our understanding of how to improve RANS modeling of this type of flowfield. 

Most importantly, transition influence was included. Two methods were tested: (1) manually zeroing out turbu- 
lence production in pre-determined regions, and (2) use of a recently published 4-equation transition model. The 
former method proved to be tedious in general, and impractical to use at high angles of attack. The latter method 
proved to be a very effective engineering tool for this complex three-dimensional high-lift flow, in spite of some con- 
vergence difficulties associated with the model. The 4-equation transition model predicted laminar regions fairly well 
in comparison with an e N method, and tended to increase lift and decrease moment, in better agreement with experi- 
ment than fully turbulent computations. Upper surface flap separation was reduced, and agreement with experimental 
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surface pressures and velocity profiles was improved. 

The current study also confirmed that including the effects of transition in combination with the brackets in the 
CFD model improved agreement with the experimental data at pre- stall angles of attack. However, the unstructured 
grids with brackets that were used had poor resolution in the wake regions above the main and flap elements. This poor 
local grid resolution yielded diffused velocity profile predictions. The precise effects of these diffuse predictions are 
not known, but it is possible that they could influence the ability of the computed flowfield to resist stall . 38 Recently- 
released experimental velocity profiles were valuable for demonstrating the influence of grid and turbulence model 
on wake shape predictions. Several different turbulence model enhancements to account for rotation and curvature 
were tested. These all generally had the desirable effects of increasing lift (reducing upper surface pressures) and 
improving the resolution of the wing tip vortex as it convected downstream. However, none of the models improved 
the predictions of the wing surface pressures near the wing tip, where the grid resolution was likely inadequate. 

Consistent and reliable prediction of the flowfield near stall angle remains elusive. There are offsetting factors 
(e.g., increasing grid resolution and accounting for transition tending to increase lift; accounting for brackets tending 
to decrease it), so that it is possible to get a reasonably good Cl, max prediction for the wrong reasons. As a result, firm 
conclusions about the efficacy of specific turbulence models near Cl, max are difficult to make without addressing all 
of the influencing factors such as transition, geometry, and grid resolution simultaneously. Accurate predictions will 
require transition modeling and inclusion of all important geometric features such as brackets, as well as adequate grid 
resolution. Particular attention needs to be paid to the wing tip and wake regions of the grid, either through a priori 
grid refinement or by using grid adaption. 
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Figure 1: Flat plate skin friction coefficient results using the 7 -Ree SST model in three different codes, M = 0.2. 
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(b) Unstructured UH16 grid (includes brackets). 


Figure 2: View of underside of NASA Trapezoidal Wing (looking forward), without and with brackets, and showing surface 
grid spacing on the main element. 
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Figure 3: View of pressure tap locations (on surface) and 7-hole-probe locations (normal to surface) over the upper surface 
of the NASA Trapezoidal Wing (pressure tap locations compared in this study are highlighted in orange). 
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(a) Lift curve, CFL3D. 


(b) Lift curve, FUN3D. 




(c) Drag curve, CFL3D. 


(d) Drag curve, FUN3D. 
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(e) Pitching moment curve, CFL3D. 
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(f) Pitching moment curve, FUN3D. 


Figure 4: Force and moment comparisons for NASA Trapezoidal Wing. 
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(a) a = 13°, flap forward station. 


(b) a = 28°, flap forward station. 
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CFL3D, SX1 (no Brkts), FT, SA 

CFL3D, SX1 (no Brkts), FT, SST 

CFL3D, SX1 (no Brkts), 7 -Re 6 SST 
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(c) a = 13°, flap at 50% span station. 


(d) a = 28°, flap at 50% span station. 
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(e) a = 13°, flap at 85% span station. 


(f) a = 28°, flap at 85% span station. 


Figure 5: Comparison of surface pressure coefficients using CFL3D on SX1 grid. 
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(a) a = 13°, flap forward station. (b) a = 28°, flap forward station. 



X 



1 1 1 1 1 1 1 1 1 1 111111 

58 60 62 64 66 68 70 

X 


(c) a = 13°, flap at 50% span station. 


(d) a = 28°, flap at 50% span station. 
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(e) a = 13°, flap at 85% span station. (f) a = 28°, flap at 85% span station. 


Figure 6: Comparison of surface pressure coefficients using FUN3D on UH16 grid. 
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(a) Main element, 15% span location, near x = 17.5 inches. (b) Main element, 83% span location, near x = 64.2 inches. 
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(c) Forward location on flap element, 83% span location, near (d) Rear location on flap element, 83% span location, near x = 

x = 69.7 inches. 75.6 inches. 


Figure 7: Comparison of velocity profiles at a = 28° using different turbulence models, on SX1/UX9 grid (no brackets). 
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(a) Main element, 15% span location, near x = 17.5 inches. 


(b) Main element, 83% span location, near x = 64.2 inches. 
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(c) Forward location on flap element, 83% span location, near (d) Rear location on flap element, 83% span location, near x = 

x = 69.7 inches. 75.6 inches. 


Figure 8: Comparison of velocity profiles at a = 28° using different turbulence models, on UH16 grid (both without and 
with brackets). 
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(a) SX1 structured grid. (b) UH16 unstructured grid. 

Figure 9: Grid section cuts at y = —73 inches (near 85% span). 




Figure 10: Residual and lift coefficient histories for a = 13° case, using CFL3D with 7 -Ree SST model. 


18 of 22 


logCy residual) 



N 


(a) Nondimensional u-velocity contours at y = —73 inches (near (b) Nondimensional u-velocity contours at y = —73 inches (near 

85% span), SST model. 85% span), 7 -Reg SST model. 
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(c) Contours of Ht / Href at y = —73 inches (color map adjusted 
to make laminar regions clearly visible), SST model. 
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(d) Contours of Ht / Href at y = —73 inches (color map adjusted 
to make laminar regions clearly visible), 7 -Reg SST model. 



(e) Upper surface streamlines, SST model. 


(f) Upper surface streamlines, 7 -Reo SST model. 


Figure 11: Comparison of SST model results with 7 -Reo SST model results, a = 13° case using CFL3D on SX1 grid. 
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(e) a = 34 °, upper surface. (f) a = 34 °, lower surface. 

Figure 12: Comparison of approximate transition location predictions; colors represent results from 7 -Ree SST model 
using CFL3D on SX1 grid (Blue=laminar, Red= turbulent); symbols represent results from Eliasson et al . 17 using e N method 
(black lines on surface represent zone boundaries). 
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Figure 13: Effects of brackets and transition on Cl at a = 13°, using FUN3D on UH16 grid. 




(a) Representative plot for configuration with no brackets. (b) Representative plot for configuration with brackets. 

Figure 14: Vorticity contours at x = 60, 70, 80, and 90 inches, a = 13°. 
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(c) Flap forward C p at a = 13 °, using CFL3D on SX1 grid. 


(d) Flap forward C p at a = 28 °, using CFL3D on SX1 grid. 



(e) Flap forward C p at a = 13°, using FUN3D (trans set) on (f) A C p between SA-RC and SA at a = 13°, using FUN3D 

UH16 grid with brackets. (trans set) on UH16 grid with brackets. 


Figure 15: Effects of various turbulence modeling corrections. 
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